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I discuss old and recent aspects of QCD jet-emission and describe how hard 
QCD results are used to construct Monte Carlo programs for generating 
hadron emission in hard collisions. I focus on the program HERWIG at LHC. 



oo : 1 The status 

^ LHC is a discovery machine, it is expected to tell us how to complete the unified theory 

O ■ of elementary interactions. New (heavy) particles are searched to indicate/confirm new 
^ . symmetries. Events with heavy particles are expected to be accompanied by an intense 
^ i emission of hadrons at short distances, and this is the domain of perturbative QCD. 
Q^l Therefore, to identify and understand non-standard events a quantitative knowledge of 
p • the characteristics of the hard radiation is strongly needed. In 1973 QCD was at the 
^ ■ frontier of particle physics (discovery of asymptotic freedom [1] and beginning of quanti- 
tative QCD studies), now in 2007 QCD is at the center of particle studies. The Monte 
Carlo programs for jet emissions [2-4] are important instruments for analyzing standard 
^ ' and non-standard short distance events. They are the Summa of most QCD theoretical 
' results and many present studies aim to improve their quantitative predictions. Thanks 
to the QCD factorization structure [5], Monte Carlo programs can be interfaced with hard 
cross sections involving also non-QCD processes (electroweak, supersymmetric, extra di- 
mension, black holes, ...). In this way, Monte Carlo generators can describe both QCD 
and non-QCD events at short distances. 

In this paper I describe the main QCD results which enter the construction of a 
Monte Carlo generator. They are so many that most of the key points will be recalled 
in a schematic way, but I hope that this short description could provide an idea of the 
reliability range of the Monte Carlo generators. For a more detailed description see [6]. 
Here, aiming to be simple and synthetic, I follow a personal point of view and the focus 
will be on the Monte Carlo event generator HERWIG [2]. Its general structure is similar 
to other important Monte Carlo generators [3,4]. In section [2] I present the scheme of 
the operations performed by Monte Carlo codes for LHC. The fact that the generation of 
events can be subdivided into successive stages is physically based on QCD factorization 
properties. The theoretical basis are discusses/recalled in section Ei In section H] I discuss 



^To appear in the volume String Theory and Fundamental Interactions, published in honour of Gabriele 
Veneziano on his 65th birthday, eds. M. Gasperini and J. Maharana, Lecture Notes in Physics, Springer, 
Berlin/Heidelberg 2007. 
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the multi-gluon soft distributions and in section [5] I describe in detail a Monte Carlo 
code for soft emissions. Although important non-soft contributions included in a realistic 
Monte Carlo are here missed, it provides a simple example containing many important 
physical effects. In section [6] I discuss non-perturbative effects which enter the Monte 
Carlo generators. The last section contains final considerations. 

2 Structure of Monte Carlo generator 

I start describing schematically the way a Monte Carlo code is organized in order to 
generate hard QCD and non-QCD events at LHC. As a specific illustration I consider the 
emission of two jets with high Et- This process is factorized into the elementary hard 
distribution, the parton densities (structure functions as in DIS) and the fragmentation 
functions (as in e~^e~): 




Elementary hard distribution 
Structure function 
Fragmentation function 



Here are the necessary factorised steps: 

• start form the hard elementary distribution Cabled with ah the incoming and cd 
the two outgoing partons. This hard distribution corresponds to QCD jet emission 
at high Et- Here one can substitute distributions for other QCD or non-QCD 
processes. There are many studies of hard distribution for processes relevant for 
LHC, see [7]. 

• generate the momenta of the hard incoming {ah) and outgoing {cd) partons (and pos- 
sible non-QCD particles). Given the hard scales Et (and possible heavy masses), the 
momenta are generated (via important sampling) in computing the total cross sec- 
tion as convolutions of the elementary distribution and the parton densities (struc- 
ture functions); 

• use the initial state space-like evolution (which at the inclusive level gives the struc- 
ture functions) to generate the "bremsstrahlung" of outgoing initial state partons 
k[,k2 - ■ ■ ■ This requires imposing a minimal transverse momentum w.r.t. the colli- 
sion direction; 

• given the outgoing hard QCD partons cd and k[,k2---, start the QCD shower 
(parton multiplication). First, from the set of these partons, identify their colour 
connections and reconstruct the set of the various primary qq dipoles. Here one 
works in the large Nc approximation so that a gluon, from the colour point of view, 
can be represented as a pair of quark-antiquark lines, a gluon is then associated to 
two dipoles; 
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• generate, for each primary dipole, the multi-parton emission according to the co- 
herent branching structure that will be illustrated in the following. This requires 
imposing a lower bound on the relative transverse momenta of final state partons 
(inside sub-jets); 

• match with the exact high order calculation, if available. It consists in weighting 
the generated event by comparing [8] the Monte Carlo distribution and the exact 
square matrix element computed to higher order [7]; 

• given the system of all emitted partons, generate the final hadrons by using a hadro- 
nisation model making hadrons out of partons. Using hadronisation models based 
on colour connections and preconfinement [9], such a process should not substan- 
tially modify [10] the structure of the hadronic radiation with respect to the partonic 
one which has been obtained in the previous steps. 

In the next sections I describe the QCD basis of these steps. 

3 The long way to Monte Carlo 

QCD has a dimensionless coupling but, even at large scale Q, when all masses can be 
neglected, the cross sections do not scale simply as powers of Q^. This is due to the 
presence of ultraviolet, collinear and infrared divergences. Ultraviolet divergences are 
responsible for the presence of the fundamental QCD scale Aqcd entering the running 
coupling. Collinear and infrared divergences are well know from QED [11]. Parton distri- 
butions can be computed only by fixing a resolution Qq (technically, a subtraction point) 
in the parton transverse momentum. Collinear and infrared divergences are responsible 
for large enhancements in these distributions which need to be resummed. Monte Carlo 
generators do actually perform these resummations as I discuss in the following. 

The possibility to resum these enhanced terms is based on specific properties of the 
collinear and infrared singularities: they factories [5,12,13]. In this way one can formu- 
late recurrence relations that lead to evolution equations. The fundamental one is the 
DGLAP evolution equation [14] resumming collinear singularities in parton densities and 
fragmentation functions. These are single-inclusive quantities, but to reach a complete 
description of an event one needs many-particle distributions so that the fully exclusive 
picture can be reconstructed (with given resolutions). The way to this is the jet-calculus 
formulated and constructed by Ken Konishi, Akira Ukawa and Gabriele Veneziano [15] 
as generalization of the DGLAP evolution equation. Therefore their work can be consid- 
ered as the basis of the Monte Carlo parton multiplication. Jet calculus leads the way 
to the evolution equation for the generating functional [12,13] of the multi-parton distri- 
butions and then to the branching probabilities for parton splitting in a way that could 
be implemented into Monte Carlo codes. The pioneering Monte Carlo codes [16-18] were 
resumming collinear singularities but only after the discovery of coherence of soft gluon 
radiation, both collinear and infrared enhanced logarithms where correctly resummed. 
The present Monte Carlo generators [2-4] fully resum not only the leading collinear and 
infrared singularities, but also relevant subleading contributions. 

In the following I describe the main theoretical points corresponding to the Monte 
Carlo steps recalled in the previous section. 



3.1 Asymptotic freedom and physical coupling 



At short distance the theory becomes free [1] and here the use of perturbation theory is 
justified. At two loops one has 

Att f 2/5i InL \ r , <5 



with (3q = 11 — |n/, Pi = 51 — yrij and nj the number of light flavours. 

To account for high order effects one needs to start from the scheme for the defini- 
tion of the running coupling. A physical definition [19] is given by the strength of the 
distribution for the emission of a soft gluon /c off a colour singlet pair of a massless quark 
and antiquark of momenta p, p. It is given by 

and corresponds to the coupling associated to the Wilson loop cusp anomalous dimension 
[20]. The relation to the MS coupling is known at three loops [21]. The argument of 
the coupling, the transverse momentum kt relative to the emitting dipole, is obtained by 
using dispersive methods [12,22] or, directly, by two loop calculations [23]. In order to 
accurately describe soft emissions, the physical coupling with the argument in ([2]) is used 
in the Monte Carlo generators. 

3.2 Coherence of soft gluons and colour connection 

Successive soft gluon emission takes place into angular ordered regions with intensities 
related to the colour charges. In the large Nc limit these regions are identified by the 
parton colour connections. To explain this one starts from the emission of a soft gluon k 
off a colour singlet qq pair, the dipole ([2]). This distribution has collinear singularities for 
6pk = or 6kp = 0. Introducing the angular variable C,ij = 1 — cosfty one can isolate the 
two singular pieces and write 



[pk)[kp) V ikp J '^'^ V ikp 

and similarly for the function '^^p{k) associated to the singularity for = 0. Performing 
the integration of over the azimuthal angle around a one has 

^%pik) = e{^pp-U), a = p,p. (4) 

This shows that the soft dipole distribution is made up of two collinear pieces, the one 
singular for k collinear to a {C,ak = 0) is (upon azimuthal averaging) bounded to a cone 
around a with opening half-angle 9pp. Since the qq dipole is a colour singlet system, the 
p and p colour lines are "connected" . 

This coherent structure can be generalized to the soft emission of a gluon A; off a 
colour singlet system made of any number of partons. Consider a, qqg colour singlet of 
momenta p and q respectively. The distribution is given by (for simplicity we take also 
the gluon q to be soft) 

Wppg{k) = Wpp{q) ■ (wpg{k) + Wgp{k) - -^Wpp{k)j . (5) 



Splitting all dipole distributions as in ([3]) one can classify all coUinear singularities in 
successive emissions within corresponding angular regions. One finds that the piece which 
is singular for k collinear to a (with a = p,p or q) is bounded to a cone around a with 
opening half-angle 6ab with b the parton colour connected to a (recall that in the planar 
limit the gluon is equivalent to a quark-antiquark pair). 

This angular ordered structure associated to colour connections at large has been 
extended [24] to the 2^2 QCD hard processes needed for LHC and used in [2]. Beyond 
large Nc, the structure of soft radiation off the 2 — > 2 hard QCD is quite more complex; 
it involves [25] rotation in the colour space for the hard matrix elements and includes 
Coulomb phase contributions. This is a very interesting contribution and would be nice 
if it could be included in a future Monte Carlo generator. 

The distribution of a soft gluon k emitted off a colour singlet pair of massive quark 
and antiquark P and P is given by 



w -(k)=-^(^-^y =^^-^^-^^ 

with (ij) = EiEj{l — ViVj cos6'jj) and Vi = ^1 — mf/Ef. While in the massless case ([3]) 
the distribution is collinear singular for k parallel to the emitting charges, in the heavy 
quark case the collinear singularities are screened: distribution vanishes for k parallel to 
the heavy quark (or antiquark) Pa and the radiation is suppressed [26, 27] in the cone 

COs6'afc > Va. 

The heavy quark screening is included into the Monte Carlo generators. One needs 
to avoid sharp cutoff around the heavy quark which, taken together with the angular 
limitations, would leave a dead cone, a phase space region without radiation. 



3.3 Sudakov form factor and jets 

An important element in Monte Carlo generator is the probability that, in a hard process, 
a parton is not radiating within a given resolution, the Sudakov form factor. To introduce 
this quantity, consider the inclusive distributions (no particle momenta are measured 
but only energy flows) which are free from collinear and infrared singularities. Classical 
examples in e"'"e~ are the jet-shapes distributions S(Q, V) with 

V = Y,v{h)- (7) 

i 

Here the sum runs over all particles in the final state (hadrons in the measurements and 
partons in the calculations). For v{k) linear in the particle momentum, such jet-shape 
observables are collinear and infrared safe. Actually individual Feynman diagrams for 
real emitted partons and virtual corrections are divergent but they are summed in such 
a way that, order by order, the infinities cancel [11] leaving finite results. 

Collinear and infrared safe jet-shape distributions Ti{Q,V) have a perturbative ex- 
pansion with finite coefficients 

S(g, V) = J^oiQ, V){1 + as(Q) c,{V) + a2(g) c,{V) + ■■■), Q » Aqcd (8) 

with T,o{Q,V) the Born distribution and Ci{V) finite functions of V expressed in terms 
of the quark, Cp, or gluon, Ca, colour charges. Actually, by inhibiting the radiation by 
taking ^ <^ 1, these coefficients are enhanced by powers of InV. A clever reshuffling 



of PT series, based on universal nature of soft and coUinear radiation (factorization) 
results [12, 13] in the exponentiated answer of the Sudakov form factor S{Q, V) 



E(Q, V) = So(g, V) ■ S{Q, V) , S{Q, V) = e'^^^'^) , 

°° (n) 
n{Q, V) = J2 {dn ln"+' V + Snln^V + ■■■). ^ ^ 

n=l 

The dn series is referred to as double logarithmic (DL) and s„ as single logarithmic (SL). 
Reliable predictions for these distributions require the matching [28] of the exact finite 
order calculation ([8]) for finite V and the Sudakov resummation iQ for small V. 

It is instructive to discuss the emergence of the powers of In V in the Sudakov form 
factor S{Q, V). They result from the incomplete cancellation of real and virtual effects. 
For V <^ 1 the real parton production is inhibited, one has v{k) < ^ 1. Since the 
virtual PT radiative contributions remain unrestricted, the divergences do cancel in the 
region v{k) < V leaving only virtual contributions for v{k) > V which produce finite 
but logarithmically enhanced leftovers. The DL contributions originate from the fact 
that each gluon emission brings in at most two logarithms (one of collinear, another of 
infrared origin). This explains the first term diln^V while the rest of the DL series is 
generated simply by the presence of the running coupling ([1]). The SL contributions, are 
necessary to set the scale of the logarithms (In" cV = In" V + n In cln"~^ V + ■ ■ ■). 

In conclusion, the Sudakov form factor S{Q, V) corresponds to the probability that 
in e'^e~ the primary quark- ant iquark pair remains without accompanying radiation up to 
resolution Qo = VQ for small V. 

To obtain the result one uses the fact that the collinear and/or infrared enhanced 
contributions factories and are resummed by linear evolution equations of the DGLAP 
type. Therefore, after factorization of collinear and infrared singularities (including soft 
gluon coherence) QCD radiation appears as produced by "independent" gluon emission 
(bremsstrahlung) . Gluon branching (into two gluons or quark-antiquark pair) enters only 
in reconstructing the running coupling ([T]) as function of transverse momentum. The 
fact that here the branching component does not contribute (within SL accuracy) can be 
understood as a result of real- virtual cancellations of singularities. Indeed, in the collinear 
limit, the transverse momentum of an emitted gluon is equal to the sum of transverse 
momenta of its decay products. Therefore, if one measures the total emitted transverse 
momentum, as in broadening for instance, it is enough to consider the contributions of 
primary bremsstrahlung gluons. Further branching does not contribute due to unitarity 
(real- virtual cancellation) . 



3.4 Structure and fragmentation functions 

Moving to less inclusive measurements one faces infinities. The simple case involves fixing 
(measuring) momentum of a hadron, e.g. that of the initial proton in DIS (structure 
function) or of a final hadron (fragmentation function), they are functions of the Bjorken 
and Feynman variables respectively 

- X, - (10) 

In DIS q is the large space-like momentum transferred from the incident lepton to the 
target nucleon P. In e^e~ annihilation q is the time-like total incoming momentum and 
P the momentum of the final observed hadron. 



In perturbative calculation, replacing the hadron with a parton, one has infinities, 
real and virtual contributions do not cancel. Soft divergences still cancel but collinear ones 
do not, making such observables not calculable at the parton level. These effects, however, 
turn out to be universal and, given a proper technical treatment, can be factored out [5] 
as non-perturbative inputs. What remains under control then is only the Q^-dependence 
(scaling violation pattern). This fact is realized in the DGLAP evolution equation which 
needs, in order to be solved, an initial condition at a low virtuality Qo- This corresponds 
to a parton resolution (or a factorized subtraction point), which absorbs all large distance 
divergences. Such "initial condition" cannot be computed by perturbative means and has 
to be provided by low scale experimental data. 

3.5 DGLAP evolution equation for DIS and e+e 

To derive the DGLAP evolution equation [14] one needs to study the phase space region 
leading to collinear singularities. The same Feynman diagrams are involved in the case of 
structure function (space- like) and fragmentation function (time- like). Therefore they can 
be studied simultaneously. First note that the Bjorken and Feynman variables ( fTOl) are 
mutually reciprocal: after the crossing operation P —P one x becomes the inverse of the 
other (although in both channels < x < 1 thus requiring the analytical continuation). 

Such a reciprocity property can be extended to the Feynman diagrams for the two 
processes and, in particular, to the contributions from mass-singularities. Consider, for 
DIS (S-case) and e~^e~ annihilation (T-case), the skeleton structure of Feynman graphs 
in axial gauge and the kinematical relation leading to the mass-singularities 




DIS or e+e skeleton graphs 



Here k[,--- k'^ are the outgoing parton systems (sub-jets). For space-like (S: < 0, ko 
entering) and time-like (T: > 0, fco outgoing) one has 

S: P^^Zi and T: Ad- ^ ^r^ . (11) 

The virtuality kf enters the denominators of the Feynman diagrams. In order for the 
transverse momentum integration produce a logarithmic enhancement, the conditions 
must be satisfied 

1^ « 1^ ^ kh « \f^\z-, (12) 

with a = —1 for DIS and a = 1 for e~^e~ . The same Feynman graphs are contributing and, 
going from S- to T-channel, the mass singularities are obtained by reciprocity: change 
z into 1/z and the momentum k from space-like to time- like. This fact is at the origin 
of the Drell-Levy-Yan relation [29] and Gribov-Lipatov [30] reciprocity which has been 
largely used in order to obtain the time-like anomalous dimensions from the space-like 
ones [31,32]. The ordering (fT2!) in the inverse fiuctuation time k'^/k+ is well known, see 
for instance [33]. 
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To make the Gribov-Lipatov reciprocity more clear, use the ordering f|T2|) in the 
computation of the probabihty D„{x^Q'^) to find a parton with longitudinal momentum 
fraction x and virtuality \k'^\ up to with cr = — 1 for the S-case and o" = 1 for the T-case. 
This ordering gives rise to the following reciprocity respecting equation [34] 

Q^dQ.D„{x, = p^z, a,) D„ g2 A ^ a = ±1 , (13) 

Jo z \z / 

with the same parton splitting kernel P{z, as) in the S- or T-channel. This equation, 
derived simply from kinematical considerations, has been (partially) tested at two [31] 
and three loop [21,34,35]. 

The reciprocity respecting equation f|T3|) is non-local since the derivative of Dcr{x, Q"^) 
in the l.h.s. involves the distribution in the r.h.s. with all virtualities larger or smaller 
than Q for cr = — 1 or (T = +1 respectively. For the use in a Monte Carlo generator 
one needs to formulate f[T5]) in terms of a local evolution equation, a Markov process. 
Formally this is easy to do: as a hard scale for the parton densities replace with 
Q\ = xQ^ in the T-case and, by reciprocity, with = x~^ in the S'-case. The 
physical meaning of these two different hard scales is well known from the studies of soft 
gluon coherence [12, 13,33,36]: in the T-case is related to the branching angle and in the 
S'-case to the transverse momentum. 

It is interesting to illustrate this. The fact that, in the T-case, the ordering variable 
is not the inverse fluctuation time /c^/fc+ (IT^ but rather the angle k'^/k'^ ~ kf/k\ ~ 6*^, 
originates from cancellations [36] due to destructive interference in the region 

T-case: z^kf < ki_^ < Ziki , (14) 

thus leaving the angular ordered region kf_^ < zfkf. Using reciprocity {zi — * 2;"^) one has 
that in the S'-case the canceling region (fT^ becomes 

S-case: \k1\ < \k'^^_^\ < zt^ | , (15) 

thus leaving the transverse momentum ordering k\^_Y < k^-. This agrees also, at small x, 
with the BFKL [37] leading order multi-parton kinematical region. 

The cancellation in the region (|T5|l has a well known physical basis for small x. 
Consider (see the skeleton graph) the successive emissions ki-2 ^i-i + and 
ki + k[ in the region <S ^ ^j-2,+ giving the leading contribution for small x. 

These cancellations result from taking into account the emission of ki off the partons fcj_2 
and k[_^ in the region f[T^ . Physically, the process can be viewed upon as an inelastic 
diffraction of the incident particle /ci_2 in the external gluon fleld of transverse size of 
order ka- In the kinematical region (fT5|) the transverse size of the parton fluctuation 
fci_2 — > ki-i + k[_i is smaller than the resolution power of the probe, /c|. In these 
circumstances the destructive interference between k^ interacting with the initial {ki-2) 
and with the flnal state (fci_i + k'-_-^) comes onto the stage. The cancellation under 
discussion is then equivalent to the general physical observation, due to V.N.Gribov, that 
inelastic diffraction vanishes in the forward direction. 

To deal with very small x one needs to resum at least all terms In" x as given 
by the BFKL equation [37] which cannot be accounted for by the collinear singularities 
resummation performed in the Monte Carlo codes. However, the evolution equation in [38] 
resums leading collinear and Inx terms (by enlarging the phase space and adding a non- 
Sudakov form factor) and allows Monte Carlo simulations [39] with the cost of generating 
events which need to be weighted. 
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4 Multi-gluon soft distributions 



CoUinear and infrared pieces of the multi-parton QCD distributions factories and can 
be reproduced by recurrence relations which can be formulated as a Markov branching 
process. This can be implemented into a Monte Carlo code and the simulation provides 
a "complete" description of the multi-parton emission in hard process. 

I illustrate in detail the case in the leading soft approximation. Although important 
non-soft contributions that are included in a realistic Monte Carlos are here neglected, 
many important physical effects are well described, in particular, large angle soft emission 
(without collinear approximation). Moreover, in this approximation the path from multi- 
gluon soft amplitudes to Monte Carlo is simple to explain. The scheme of the presentation 
involves the following steps: 

• multi-gluon soft distributions. They are computed in the leading soft approximation 
and in the planar approximation; 

• recurrence relation for the multi-gluon soft distributions. This is obtained by intro- 
ducing the generating functional for all multi-gluon distributions [12] and deriving 
the evolution equation. From the generating functional one computes observables 
as it will be discussed in subsection 14. 2[ For collinear and infrared safe observables 
such as jet-shape distributions the cutoff contributes only with power corrections; 

• Markov process and Monte Carlo implementation. Here one needs to include proper 
cutoff for collinear and infrared singularities. This will be discussed in the next 
section [51 

• from parton to hadron emission. This will be discussed in section O 

The starting point is the amplitude for the emission of n soft gluons gi, ■ ■ ■ ,qn off 
a primary colour singlet qq pair of momentum p, p. It is represented as a sum of Chan- 
Paton factors with the coefficients given by colour- ordered amplitudes. We consider the 
contribution with a single Chan-Paton factor (topological expansion [40]) 

Mnippqi ■■■qn) = J^IA'^'i ■ ■ ■ A'^- M^ipq,, ■ ■ ■ q,^p) , (16) 

the sum is over the permutation 7r„ of colour indices. A" are the SU {N^) matrices in the 
fundamental representation. The softest emitted gluon q^ factorizes and one has [12,42] 

M„(... £m£'- ■■)= ^7sM„_i(- ■■£/■■•)■ f 7^-^^) . (17) 

\{mm) {qe'qmjj 

The softest gluon is emitted by the two partons neighbouring in colour space. This 
approximation is accurate in the soft limit without any collinear approximation. From this 
factorized structure one deduces a recurrence relation and computes all colour-amplitudes 
in the soft limit. Summing over the polarization indices, the squared averaged colour- 
amplitude is given, for the fundamental colour permutation, by 

\Mr,ipq,---qnp)\'=\Mo\'{2g',rWpp{q^---qn), ^^^^^^^ " " " = (^^J^T^^ • ^^^^ 

This very simple result for the square amplitude is valid for any energy ordering and 
depends only on the colour ordering. Note that here one takes the square of the same 



colour-ordered amplitude. Indeed M„(7r^)M*(7r„) with 7r„ and tt^ two different colour 
permutations cannot be expressed in a closed form for any n. On the other hand contribu- 
tions from different permutations enter the calculation of the averaged squared amplitude 
|A1„p. A close expression for this distribution for any n is obtained only in the planar 
approximation [41]. To see this observe that 



Tr(A.„A.r) = 2C^^^J V^"iV^J ^^^^ 

with Aff„ = {A^^-'-A""} and A^t = {A°" ■ ■ ■ A"^}. Taking instead two different colour 
permutations one has that Tr(A7r; \t,t) is suppressed at least by l/N^- Therefore, only in 
the planar approximation one can use the simple result in (fT8|) and obtains [12] 



\Mn\' = ^ (N^gir E ^pp(^^^ ■ ■ ■ (20) 



where ctq = 2Cf\Mq\'^ and symmetrisation has been taken into account. 

The distributions ( |T8i) contain the leading infrared singularities: for any colour per- 
mutation one has Wpp ~ (cui ■ ■ ■ujn)"'^ with Ui the energy of gluon g^. They contain also 
the leading collinear singularities for 6ij = with ij two partons neighbouring in colour 
(thus there are up to n colhnear singularities). 

An alternative way to obtain the the multi-gluon colour amplitude is based on the 
helicity techniques [43]. For qq with + and — polarization, the leading soft contribution 
is obtained when all gluons have + helicities and the recurrence relation (fTTI) reads (for 
opposite helicities the result is the complex conjugate one) 

M„(- ■ ■ £m f ■ ■ ■ ) = M„_i(- ..&'...). ^^^f^ , {qq') = ■ e^V , (21) 

{qiqm) {qm.qe) 

with qm the softest gluon, z the longitudinal direction and the phase 



e^^=y^^--^j, <l. = q. + ^qy. (22) 

The solution of this recurrence for the amplitude is very simple; it is the same for any en- 
ergy ordering and depends only on the colour ordering. For the fundamental permutation 
one has 

Mn{pqi ■ ■ ■ qnp) = g: Mo , . , (23) 

{pqi} ■ ■ ■ \qnp) 

with squared amplitude given by (fT8l) . This shows the well known result that non-planar 
contributions, obtained from M„(7r„) ■ M'^{ti'^ for two different colour orderings, have the 
same soft singularities but reduced number of collinear singularities. 



4.1 Virtual correction, generating functional and evolution 

To compute observables one needs to supplement the multi-gluon soft distributions fl20l) 
with the related virtual corrections. For infrared and collinear safe observables, such as 
jet-shape distributions, the infrared and collinear singularities in 0181) has to be canceled 
by corresponding singularities in virtual corrections. One way to compute the virtual 
corrections, at the same level of accuracy in the soft limit as for real emission contribution. 



in 



consists of performing the integration over the virtual gluon energy by the Cauchy method 
and then taking the soft hmit for the virtual gluon. This way one also regularizes the 
ultraviolet divergences by neglecting the divergent contribution from the contour at the 
infinity of the complex energy plane. By properly choosing a constant this regularization 
corresponds to the physical scheme in ([1]). The virtual corrections so computed can be 
included into the generating functional for the multi-gluon soft distributions. The result of 
this study not only gives the relevant virtual corrections but, due to the simple structure 
of (l20l) in the planar approximation, gives the branching structure of multi-gluon soft 
emission leading to the Monte Carlo generator. 

Consider the soft distribution da^^^ for the emission of n gluons off a colour singlet 
dipole ab (thus one generalizes the primary dipole pp to a general dipole with a and b 
in arbitrary directions). For each emitted soft gluon one introduces a source function 
u{qi) and defines the generating functional as 

Gab[E,u] = J2^ (24) 

with E = Q/2 the hard scale. This functional depends on the directions a and b of the 
primary dipole. By setting all u{qi) = 1 one has Gab[E,l] = 1. Using (I2U|) one has the 
real emission contribution for the generating functional 

G^f [E, u] = J2 [ n("^^(^*) ^^d^^QiE-iU,)\ ■ Wabiqi ' ' ' qn) , (25) 

with Us = Ncas/ir. Here one neglects corrections (planar limit) and uses the soft 

approximation for the phase space Ui <^ E. Symmetry of the phase space is used. The 
condition Gab[E, 1] = 1 must be satisfied only after including the virtual corrections. To 
include them we construct the evolution equation for the generating functional. To this 
end we use the fact that the very simple expression flTSl) has the following factorization 
property 

Wab{qi ■■■qn) = Wabiqe) ■ Waiiqi ■ ■ ■ qe-i) ■ Wib{qm ■■■qn) , (26) 

with qe one of the soft gluons and Wab{q) the dipole distribution ([3]). Taking qe as the 
hardest (soft) gluon and differentiating fl25l) with respect to E, thus setting Ui = E, one 
obtains [44] 

EdEGab[E,u] = J ^^^^^u{q)Ga,[E,u]-G,b[E,u]-G,b[E,u]], (27) 

with = 1 — cos 6ij. The negative term in the integrand originates from the virtual 
corrections obtained via Cauchy integration as mentioned before. Since they are evaluated 
within the same soft approximation used for the real contributions, at the inclusive level 
they cancel against the real contributions giving the correct constraint Gab[E, 1] = 1- Both 
the real emission (first term in the integrand) and the virtual correction (second term) are 
coUinear and infrared singular. For inclusive observables, (i.e. for suitable sources u{q)) 
these singularities cancel. This evolution equation accounts for coherence of soft gluon 
radiation [12, 13]. 

4.2 Observables in the soft limit 

Using Gab[E,u] one obtains all inclusive distributions in the soft limit. No coUinear 
approximations are involved in (!20!) . therefore the functional Gab[E,u] gives quantities 
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which involves also large angle soft emission. Let me first recall some observables which 
are coUinear singular around the primary partons a and h. 



Collinear observables. The simplest one is the multiplicity of soft gluons with reso- 
lution Qq. Taking u{q) = u this observable is defined as, see ([2] 



(n) 

nab{E) = duGab{E, u) = V . (28) 



It is easy to derive from fl27|) the well known result [36] for the multiplicity 



ME).„i?exp||^^^ (29) 

with n^^fj the non-perturbative initial condition. Similarly one derives the fragmentation 
function Dab{x, E) by taking the source u{q) = u{x) with x the soft gluon energy fraction 

Dab{x,E) = -^Gab[E,u] . (30) 

dU[X) u(x)=l 

Soft gluon coherence here is shown by a depletion of radiation [12, 13] at small x. 

Observables at large angle. The simplest case is the distribution discussed in [45] 
of heavy systems of mass Ai emitted in e^e~ at large angle p = ^(1 — cos^) and small 
velocity. The heavy system (typically a heavy qq system) originates from a gluon in the 
cascade. The collinear singularities are screened by Ai so this distribution is finite and 
given by a function of the SL quantity 

dqt 

r = — as{qt) ■ (31) 
Jm 

It is interesting that this distribution I{p, r) satisfies an equation with a structure similar 
to the BFKL equation [45] and then its asymptotic behaviour in r involves the BFKL 
characteristic function. One has 

nP.r)~'^-'^e''4g^. fl = 28C(3). (32) 

The functional Gab[E,u] is suited to give the distributions in the energy emitted away 
from jets. Such distributions do not have collinear singularities, but only infrared ones. 
An example in e~^e~ is the distribution in energy recorded outside a cone around the 
thrust (this is a typical "non-global" jet observable [46]): 

E(E, E„,,) = x: / e ( - y: ku 

Since the jet region is excluded, there are no collinear singularities to SL accuracy and 
the resummed PT contributions come from large angle soft emission. Here resummation 
is complex but informative. It brings information on the QCD radiation between jets, a 
region interesting for understanding colour neutralization among jets. 





out 




thrust 1 Oin\ 






axis \ . 1 








out 
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It is interesting to discuss this quantity in some detail since it illustrates the structure 
of ( 1271) . First observe that the distribution depends on E and Eout through the SL function 
r given by ( l3Ti) with M. — -Eout- To obtain S(r) from Gab[E,u] one takes u{q) = 
away from jets and u{q) = 1 inside the jet region. From (1271) one derives the evolution 
equation [44] 

drT.ab{r) = -Sab^ab{r) + [ ^ ^^^^ | S^g (t) ■ Sgfe (t) - ^^^(r) ) , (33) 

Jin 47r ^aqt.qb ^ > 

with Sah related to the Sudakov form factor 

S{r) = e-^^- , Sab 



[ (34) 



Equation (I55|) has a bremsstrahlung (first) and branching (second term) components: 






hremsstrahlung component branching component 

The hremsstrahlung component resums contributions from gluons emitted in the recorded 
region outside the cone. These contributions are the only ones present for the global jet 
observables considered in the previous subsection. Here, since the coUinear singularities 
are screened by the cone ^in, the Sudakov form factor is a SL function. 

The branching component resums contributions from gluons emitted inside the jet 
region. These gluons need to branch in order to generate decay products entering the 
recorded region. Here real-virtual cancellation is incomplete and virtual enhanced con- 
tributions are dominating thus leading to a strong suppression of the distribution which 
asymptotically turns out to be Gaussian in r. 

The Monte Carlo generator [2] resums only collinear singularities therefore it does not 
fully resum soft emissions at large angles although phenomenologically, it turns out [47] 
that the most important pieces are correctly reproduced due to soft gluon coherence. 



5 Monte Carlo simulation for soft emission 

The evolution equation (|27|) can be formulated as a Markov process and then numerically 
solved. This Monte Carlo procedure has been introduced in [46] to study non-global 
distributions. A similar procedure based on dipole branching is used in the Monte Carlo 
generator [4]. 

To construct a Monte Carlo generator from (1271) one splits the real and virtual cor- 
rections. To do so it is necessary to introduce a cutoff Qq in transverse momentum (the 
argument of Us) giving the Sudakov form factor 

lnSa,{E) = -r^ ■eig.ab-Qo), cil, = 2ul^-^ (35) 

JQo J Uq^qb kab 

which is the solution of (i27l) with the real emission piece neglected. Here qtab is the 
transverse momentum of q with respect to the a6-dipole. Then the evolution equation 
( 127|) can be integrated to give (the cutoff Qo dependence is implicit) 

Gab[E] = SabiE, Qo) + / dVab{E, Ug, Qg) U{q) Gaq[uJg, U] ■ Ggb[ujg, u] , (36) 



where one has introduced the probabihty for dipole branching: [ab) [aq) {qb) 
,^ o^ ( dujq SabiE)\ ( dilg asUb\ ^^ 

dVab{E,ujg,Qg) = \ — ^ . . > < -r^ -r-r r ' ^w^^-^o) • (37) 

UJq Sab{uJg) J t 4:71 Uq^qb ) 

To see how this could be used in a Monte Carlo simulation one writes dVabiE,uj,Q) in 
the equivalent form (the bound qtab>Qo is implicit) 

dVab{E, u, n) = drab{E, uj) ■ dRab{n) (38) 

with 

rab{E,ujg) 

dRabi^q) 

The integral of the branching probability gives 

J dVab{E,UJ,n) = l-Sab{E), (40) 

and this shows that the Sudakov factor Sab{E) gives the probability for not emitting a 
gluon within the resolution Qq in qtab- 

The probability distribution dVab{E,uj,fl) can be used to generate Monte Carlo 
events distributed according to QCD in the soft and planar approximation. Using sets of 
random numbers < p < 1 the procedure is the following: 

1. take the afe-dipole with the energy scale E and compare the Sudakov factor Sab{E) 
with p. If p < Sab{E) then the a6-dipole does not emit any soft gluon within the 
resolution. In the opposite case the dipole is emitting a soft gluon with energy ujg 
given by solving the equation p = rab{E,uJq); 

2. obtain the direction Qq by sampling the distribution dRab{Qq)- At this point, from 
the a6-dipole one has generated two dipoles: aq and qb, both at the new energy 
scale LOq] 

3. repeat the procedure for each new generated dipole till no dipole emits any more 
within the resolution. 

At the end of this procedure one is left with a Monte Carlo event: a collection of 
emitted soft gluons qi - ■ ■ qn together with the primary partons a, b. These events are 
distributed with the QCD probability so they ca be used to compute any soft distribution 
as discussed in subsection 14.21 

Such a Monte Carlo simulation, based on evolution equation in energy, is then a 
successive emission of softer and softer gluons. Angles are given by the dipole distribu- 
tion ([3]) so they are ordered (upon azimuthal average) and coherence is automatically 
implemented. 

In order to obtain a realistic simulation one needs to overcome the soft approximation, 
that is, to take into account the recoil in the emission and the non-soft pieces of the gluon 
splitting function (only the singular pieces are present in fl271) ) 

Pg^ggiz) = nJ- + -^ + zil-z)-2] . (41) 



SabjE) 
Sabi^q) ' 



ab 



ab 



dVtq ttsi 
47r ^aq^qb 



drab{E,UJq) = 1 - Sab{E) 



dRab{Qq) = 1 



(39) 
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Similarly, one needs to account also for the quark branching channels. All these points 
are accounted in the present realistic Monte Carlo generators. Their basis is an evolution 
equation in angle rather than in energy (as ( 1271) ). However this implies that one considers 
collinear approximations in the emission thus soft radiation at large angles are not fully 
accounted for. 

6 From partons to hadrons 

The above description of the Monte Carlo code refers to the generation of events with 
emission of partons (possibly together with non-QCD particles) which, due to the presence 
of collinear and infrared singularities, requires a cutoff Qo- The main questions are then: 
how to go from partons to hadrons and how much a phenomenological hadronisation 
model affects and distorts the QCD radiation generated perturbatively. A suggestion 
on hadronisation models which do not substantially modify the peturbative radiation is 
provided by preconfinement [9]. 

Preconfinement. The basis is again the Sudakov function which suppress the proba- 
bility of "non-emitting". Consider, in the planar approximation, two colour connected 
partons emitted in a hard collision at scale Q and with resolution Qq. Colour connection 
means that the quark colour line of one parton ends into the antiquark colour line of 
the other parton (in the planar approximation a gluon could be, from the colour point 
of view, described as a pair of qq colour lines). Thus no gluons are emitted within the 
resolution Qq by this colour line and a Sudakov form factor arises which forces the two 
colour connected partons to form a system of mass of order Qq (even for very large Q). 
The system of the quark and antiquark in question forms a colour singlet of small mass. 
Although this is not yet an indication of confinement (the colour system should be local- 
ized in space), such a preconfinement property suggests that any hadronisation models 
that associates hadrons to colour connected partons would not distort the perturbative 
structure of the QCD radiation: parton and hadron flows are similar within the resolu- 
tion Qq. Preconfinement is then related to the property of local hadron-parton duality [10] 
which has been phenomenologically well tested. 

Power corrections. Other non-perturbative effects are the power corrections to the 
observables. They result from the non-convergence the PT expansions even if the co- 
efficients are finite as in ([8]) and (Q. As a consequence all PT predictions are affected 
by corrections in powers of Aqcb/Q with coefficients determined by NP effects. An im- 
portant NP effect, present in short distance quantities, is that the running coupling is 
involved at any scale smaller than Q. For example, the average value of in ([7]) is given 
by an integral of the type 



where the virtual momentum kt in the Feynman diagrams runs into the large distance 
region (although the observable is dominated by short distance physics). Since the ob- 
servable is collinear and infrared finite, for kt ^ the Feynman integrand is regular 
{y{kt/Q) ~ kt/Q) so that the integral is finite, apart from the presence of as{kt) which 
enters the confinement region. Mathematically this is reflected into the fact that, although 




(42) 



all PT coefficients in as{Q) are finite, the expansion is non- convergent [48] (renormalon 
singularity) . 

The fact that the running coupling entering the NP region is at the origin of the 
leading power correction can be checked phenomenologically. From the study of jet-shape 
observable one finds [49] that, within 10 — 20%, the power corrections are described by 
the same parameter accounting for the running coupling in the NP region. In the Monte 
Carlo generators one sets a cutoff Qo in the argument of the coupling and this does bring 
in these physically relevant power corrections at the perturbative — parton — stage. 
Instead, power behaving contributions to jet shapes arise at the hadronisation level [51]. 

Underlying event. Another important NP component in the Monte Carlo for LHC 
is the presence of radiation besides the one emitted in the hard event. This is typically 
around the beams as for the peripheral interactions (events at low Et)- Perturbative 
QCD does not provide indication for this component. Thus there are various models 
which needs to be studied [50] at the Tevatron together with the extrapolation at LHC. 

7 Conclusion 

What I have discussed shows that the Monte Carlo generators involves the entire Summa of 
hard QCD results and provide a framework for many future QCD and non-QCD studies. 
The general attempts to improve the Monte Carlo generators go in the directions of 
making the quantitative predictions both more reliable (by adding new theoretical QCD 
results and phenomenological studies) and more general (by including also electroweak 
and beyond the standard model physics). As far as the ffist direction, I have mentioned 
the works made to include in the Monte Carlo generator the known exact higher order 
distributions [8] . As also mentioned, it is interesting to include into the present generators 
reliable predictions on large angle soft emission (see subsection 14. 2p . This would require 
also the need to account for non-planar corrections by studying colour rotations involved 
in the colour structure of ensembles of more than three hard partons (see [25]). 

The three key elements in a Monte Carlo generator for jet emissions are the QCD 
factorization properties, the branching algorithm and the procedure for converting par- 
tons into hadrons. As I have mentioned, Gabriele Veneziano has either contributed to or 
started each of these three key developments: The Monte Carlo generators are based on 
factorization of QCD collinear singularities [5]. Jet calculus [15] leads to the evolution 
equation for the generating functional for multi-parton distributions which can be formu- 
lated as a Markov process. Moreover, the preconfinement property [9] is at the basis of 
hadronisation models that do not destroy the QCD radiation structure. 
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